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OBJECTIVES 


1.1  Original  Objectives 

It  was  proposed  to  conduct  a  three-year  research  effort  focused  on  developing,  validating,  and 
exploiting  a  high-fidelity  and  high-performance  numerical  simulation  technology  for  predicting  the 
effects  of  combat  damage  on  the  flight  and  flutter  envelopes  of  military  aircraft.  To  this  effect,  it 
was  proposed  to  build  on  previous  achievements  in  CFD  (Computational  Fluid  Dynamics)-based 
nonlinear  computational  aeroelasticity  funded  by  AFOSR  (Grant  F49620-99-1-0007),  and  perform 
the  following  research  tasks: 

•  Advanced  modeling  and  computational  methodologies.  It  was  proposed  to  research  and  develop 
high-fidelity  and  high-performance  methodologies  for  modeling  and  simulating  the  effect  of 
damage  on  the  aerodynamic  performance,  stability,  control,  and  aeroelastic  behavior  of  mili¬ 
tary  aircraft.  These  include,  among  others,  the  concept  of  “phantom”  elements  for  facilitating 
the  introduction  of  damage  in  the  computational  fluid  and  structure  models,  fast  methods 
for  evaluating  transient  effects,  and  the  sensitivity  analysis  by  the  adjoint  method  of  flight 
mechanics  and  flutter  criteria  such  as  lift,  drag,  stability  derivatives,  and  flutter  speed  with 
respect  to  various  damage  parameters  including  the  size  and  location  of  a  hole. 

•  Dynamic  live-fire  ground-testing  methodology.  Assisted  by  the  developed  simulation  technol¬ 
ogy,  it  was  also  proposed  to  investigate  a  reconfigurable  wing  loading  methodology  based 
either  on  tethers  or  embedded  piezoelectric  actuators  which  replicates  in-flight  loading  con¬ 
ditions,  and  allows  the  structure  to  properly  react  when  damaged. 

•  Validation.  Finally,  it  was  proposed  to  collaborate  with  Dr.  Greg  Czarnecki  and  his  team  at 
the  Wright-Patterson  Air  Force  Base  (AFB)  to  validate  the  proposed  simulation  technology, 
and  assess  the  feasibility  of  the  proposed  dynamic  loading  methodology.  It  was  also  proposed 
to  work  with  Dr.  Mike  Love  at  Lockheed-Martin  Aeronautics  to  explore  how  the  proposed 
modeling  and  simulation  technology  can  help  in  reducing  some  of  the  test  costs  of  the  JSF 
(Joint  Strike  Fighter)  program. 

1.2  Revised  Objectives 

After  one  year  of  performance  on  this  grant,  the  Flight  Test  Center  at  the  Edwards  Air  Force  Base 
requested  an  additional  emphasis  on  aeroelastic  reduced-order  models  to  speed-up  the  prediction  of 
the  flutter  envelopes  of  modern  aircraft  fighters.  This  request  was  accommodated  in  the  form  of 
an  additional  task  (see  Research  Task  8  below). 

2  STATUS  OF  EFFORT  AND  ACCOMPLISHMENTS 

During  the  third  year  of  funding,  effort  focused  on  Research  Task  6  and  Research  ffhsk  8. 
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2.1  Outcome  of  Research  Task  6:  A  Stress-Control-Based  Live-Fire  Ground 
Testing  Methodology 

LFT  is  supposed  to  be  conducted  as  if  the  aircraft  was  in  flight,  and  had  been  hit  by  an  anti-aircraft 
artillery  round.  Typically,  it  is  performed  on  a  wing.  The  tanks  are  loaded  with  fuel,  and  high- 
velocity  air  is  blown  by  a  battery  of  jets  across  the  wing  which  is  loaded  by  computer-controlled 
hydraulic  jacks  to  simulate  in-flight  loads.  In  one  test,  an  explosive  bullet  is  fired  away  from  the 
fuel  tank.  In  another  test,  it  is  fired  inside  the  fuel  tank.  The  explosive  bullet  generates  a  shock 
wave  that  travels  through  the  fuel  and  imparts  loads  on  the  wing’s  skin  and  internal  structure. 
Consequently,  a  portion  of  the  wing  skin  deforms  into  the  air  stream  then  rips  off.  If  the  damaged 
wing  remains  largely  intact,  it  is  currently  concluded  that  the  limited  structural  damage  would 
enable  a  pilot  to  fly  the  airplane  home  [1]. 

The  LFT  procedure  outlined  above  has  at  least  one  main  deficiency;  the  high-velocity  air  blown 
by  the  jets  and  the  computer-controlled  hydraulic  jacks  applied  to  the  tip  of  the  wing  do  not 
reproduce  the  in-flight  loads.  In  particular,  the  hydraulic  jacks  introduce  an  unrealistic  restraint 
at  the  tip  of  the  wing  whose  root  is  already  cantilevered. 

To  reproduce  the  in-flight  loads,  a  loading  methodology  should  reproduce  the  true  stress  state 
of  the  structure.  Given  a  set  of  flight  conditions,  this  stress  state  can  be  predicted  fairly  accurately 
by  numerical  aeroelastic  simulation.  Hence,  the  main  idea  reported  here  is  to  consider  a  wing, 
instrument  it  with  a  number  of  embedded  actuators,  and  program  these  actuators  to  produce 
the  desired  aeroelastic  stress  state  field  denoted  here  by  s.  Such  actuators  could  be,  in  principle, 
piezoelectric  actuators,  thermal  actuators,  simple  tethers,  or  any  combination  of  these.  Unlike 
hydraulic  jacks,  embedded  actuators  do  not  restrain  the  wing  at  any  location,  and  potential  tethers 
soften  the  restraints  they  induce  because  of  their  compliance. 

Let  gi  denote  the  gain  of  the  i-th  embedded  actuator,  and  Uj  denote  the  displacement  field  of 
a  wing  structure  due  to  this  actuator  and  its  gain  gi.  Under  the  usual  linear  assumption,  when  all 
Na  embedded  actuators  are  activated,  the  total  displacement  of  the  wing  is 

i=Na 

i=I 

where  U  is  the  matrix  of  displacement  fields  Uj,  and  g  is  the  vector  of  gains  gi.  From  Eq.  (1), 
it  follows  that  if  s,  denotes  the  stress  field  associated  with  Uj,  then  the  stress  state  of  the  wing 
structure  is  given  by 

i=Na 

s  =  ^  Sigi  =  Sg 

i=l 

where  S  is  the  matrix  of  stress  fields  s,  associated  with  the  gains  gi. 

Hence,  the  crux  of  the  idea  reported  here  is  to  find  the  gain  vector  for  which  s  =  s  —  that  is, 

Sg  =  s.  (2) 

Unfortunately,  Eq.  (2)  can  be  exactly  satisfied  only  if  the  total  number  of  actuators  Na  is  equal 
to  the  number  of  degrees  of  freedom  of  interest,  and  using  such  a  large  number  of  actuators  is 
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unfeasible  as  well  as  undesirable.  In  particular,  the  component  of  the  structure  that  is  expected  to 
be  damaged  should  not  be  instrumented  so  that  no  actuator  is  destroyed  and  the  in-flight  loads  are 
reproduced  as  much  as  possible  after  damage  is  inflicted  on  the  instrumented  structure.  Hence,  a 
residual 

r  =  s  -  Sg  (3) 

is  to  be  accepted  but  minimized.  Here,  the  Euclidean  metric  is  chosen  for  this  purpose.  Further¬ 
more,  it  is  well-known  that  some  actuator  locations  are  more  efficient  than  others  in  prescribing 
a  particular  shape  [2,3],  and  therefore  a  particular  stress  state.  For  this  reason,  a  good  strategy 
consists  in  minimizing  |  |s  —  Sg|  [2  over  both  the  locations  and  gains  of  the  Na  actuators,  for  different 
values  of  Na- 

Alternatively,  after  damage  is  inflicted  on  the  wing,  it  can  be  inspected  and  described  by  holes 
in  the  computational  structural  dynamics  and  computational  fluid  dynamics  models.  Then,  the 
vector  of  actuator  gains  can  be  reset  to  the  value  which  minimizes 

||rl|2-||sd-Sg||2,  (4) 

where  is  the  stress  state  predicted  for  the  damaged  wing  by  numerical  aeroelastic  simulation. 
Once  the  damaged  wing  is  loaded  according  to  this  reconfigured  gain  vector,  its  structural  integrity 
can  be  assessed,  and  further  numerical  simulations  can  be  performed  to  evaluate  its  remaining  flight 
capability. 

In  summary,  the  stress-control-based  live-fire  ground  testing  methodology  developed  so-far  is 
governed  by  the  mathematical  problem 


ggRWo 

(5) 

C(g)  <  0, 

(6) 

where  C(g)  is  a  matrix  of  constraints  specifying,  for  example,  that  each  actuated  member  does 
not  exceed  a  certain  percentage  of  the  yield  stress  and  that  the  total  energy  of  actuation  does  not 
exceed  a  certain  threshold.  Note  however  that  the  solution  of  the  above  minimization  problem  is 
adequate  for  the  envisioned  application  if  and  only  if  s  G  .R(S),  and  that  the  range  of  S,  R{S), 
depends  on  the  number,  location,  and  type  of  the  actuators. 

Let  Sq  denote  the  matrix  of  stress  fields  Sj  associated  with  the  gains  gi  when  all  members  of 
the  test  structure  are  actuated.  If  s  ^  H(Sa),  this  “range  check”  indicates  that,  using  the  chosen 
type  of  actuators,  the  test  article  cannot  be  put  in  the  aeroelastic  stress  state  §  even  if  each  one  of 
its  members  is  instrumented.  In  this  case,  the  range  check  also  suggests  that  additional  internal 
members  and/or  external  members  such  as  tethers  are  needed  to  control  the  stress  sate  of  the  test 
article.  Then,  if  additional  members  are  introduced  in  the  structure  for  controlling  its  stress  state, 
stilfness-and-mass  checks  must  then  be  performed  a  posteriori  to  ensure  that  the  few  additional 
members  have  not  significantly  changed  the  characteristics  of  the  structure  —  for  example,  its  eigen 
modes. 

Advances  in  the  methodology  outlined  above  including  the  use  of  one-dimensional  versus  bender 
actuators,  its  accuracy,  feasibility  (no  member  yielding,  minimum  and  practical  energy  require¬ 
ment),  and  demonstration  for  the  ARW2  wing  have  been  reported  in  2004  in  the  AIAA  paper  [7] 
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and  documented  in  the  attached  comprehensive  Ph.  D.  report  entitled  “A  Stress-Control-Based 
Live-Fire  Ground  Testing  Essentially,  it  was  found  that  at  least  for  the  ARW2  wing,  tethers  are 
almost  always  required  to  reproduce  the  in-flight  loads  with  sufficient  accuracy. 

In  order  to  check  whether  the  above  result  is  primarily  due  to  the  fact  that  the  ARW2  wing  has 
a  large  aspect  ratio,  the  proposed  loading  methodology  was  next  assessed  for  de/ta-like  wings.  The 
smaller  aspect  ratio  of  such  wings  tends  to  make  their  structures  statically  indeterminate,  which 
in  turn  tends  to  improve  the  performance  of  self-equilibrated  one-dimensional  actuators.  More 
specifically,  two  such  wings  were  considered:  a  high-speed  civil  transport  wing  and  an  F-16  wing  in 
clean  configuration.  In  both  cases  and  for  several  damage  scenarios,  it  was  also  found  that  external 
tethers  cannot  be  avoided  if  in-flight  loads  are  to  be  reproduced  accurately.  More  specifically,  it  was 
found  that  the  entire  objective  of  the  loading  methodology  can  be  achieved  using  only  a  few  tethers. 
For  example,  it  was  found  that  for  the  F-16  wing,  eight  well-positioned  tethers  can  reproduce  in 
most  cases  the  in-flight  stress  states  of  the  wing  with  a  global  relative  error  of  the  order  of  10%. 
On  the  other  hand,  it  was  also  found  that  the  current  usage  of  computer-controlled  hydraulic  jacks 
at  the  tip  of  a  wing  to  simulate  in-flight  loads  typically  results  in  relative  errors  that  exceed  the 
50%  level. 

2.2  Outcome  of  Research  Task  8:  Reduced-Order  Modeling  for  Near- Real-Time 
Flutter 

In  order  to  enable  real-time  flutter  analysis,  a  computational  methodology  based  on  adaptive 
reduced-order  modeling  (ROM)  was  developed.  In  this  methodology,  the  structure  is  represented 
by  a  truncated  modal  set,  the  fluid  by  a  reduced-order  basis  obtained  by  Proper  Orthogonal  De¬ 
composition  (POD)  for  a  given  Mach  number,  and  the  resulting  basis  is  adapted  for  variations 
in  the  Mach  number  by  interpolation  of  the  subspace  angles.  This  computational  technology  was 
implemented  in  AERO  and  applied  to  a  complete  F-16  fighter  configuration.  Good  correlations 
between  full-order  simulation  results,  reduced-order  simulation  results,  and  flight  test  data,  were 
obtained  and  documented  in  an  AIAA  paper  that  appeared  in  2005. 
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5.2  Transitions 

The  second-generation  nonlinear  aeroelastic  simulation  platform  AERO  developed  by  the  Principal 
Investigator  and  his  research  group  was  transitioned  to  the  Flight  Test  Center  at  the  Edwards  Air 
Force  Base  (Point  of  Contact/User:  Dr.  John  Sun)  for  application  to  the  aeroelastic  modeling 
and  simulation  of  the  JSF.  It  was  also  delivered  to  the  Naval  Research  Laboratories  (Point  of 
Contact/User;  Dr.  John  Michopoulos)  for  enabling  research  in  the  design  of  a  dynamic  data-driven 
system  for  structural  health  monitoring  and  critical  event  prediction. 
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